Electronic transport properties of graphene 
nanoribbons 



Katsunori Wakabayashi^'^, Yositake Takane^, Masayuki 
Yamamoto^ and Manfred Sigrist^ 

^International Center for Materials Nanoarcliitectonics(MANA), National Institute 
for Materials Science (NIMS), Namiki 1-1, Tsukuba 305-0044, Japan, 
^PRESTO, Japan Science and Technology Agency(JST), Kawaguchi 332-0012, Japan, 
■^Department of Quantum Matter, AdSM, Hiroshima University, Higashi-Hiroshima 
739-8530, Japan, 

^Theoretische Physik, ETH-Ziirich, Ziirich CH-8093, Switzerland. 

Abstract. We will present brief overview on the electronic and transport properties 
of graphene nanoribbons focusing on the effect of edge shapes and impurity scattering. 
The low-energy electronic states of graphene have two non-equivalent massless Dirac 
spectrum. The relative distance between these two Dirac points in the momentum 
space and edge states due to the existence of the zigzag type graphene edges are decisive 
to the electronic and transport properties of graphene nanoribbons. In graphene 
nanoribbons with zigzag edges (zigzag nanoribbons), two valleys related to each Dirac 
spectrum are well separated in momentum space. The propagating modes in each 
valley contain a single chiral mode originating from a partially flat band at band center. 
This feature gives rise to a perfectly conducting channel in the disordered system, if 
the impurity scattering does not connect the two valleys, i.e. for long-range impurity 
potentials. Ribbons with short-range impurity potentials, however, through inter- 
valley scattering display ordinary localization behavior. On the other hand, the low- 
energy spectrum of graphene nanoribbons with armchair edges (armchair nanoribbons) 
is described as the superposition of two non-equivalent Dirac points of graphene. In 
spite of the lack of well-separated two valley structures, the single-channel transport 
subjected to long-ranged impurities is nearly perfectly conducting, where the backward 
scattering matrix elements in the lowest order vanish as a manifestation of internal 
phase structures of the wavefunction. For multi-channel energy regime, however, 
the conventional exponential decay of the averaged conductance occurs. Symmetry 
considerations lead to the classification of disordered zigzag ribbons into the unitary 
class for long-range impurities, and the orthogonal class for short-range impurities. 
Since the inter- valley scattering is not completely absent, armchair nanoribbons can 
be classified into orthogonal universality class irrespective of the range of impurities. 
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1. Introduction 

Recently graphene, a single-layer hexagonal lattice of carbon atoms, has emerged as a 
fascinating system for fundamental studies in condensed matter physics, as well as the 
promising candidate material for future application in nanoelectronics and molecular 
devices. [1] The honeycomb crystal structure of single layer graphene consists of two 
nonequivalent sublattices and results in a unique band structure for the itinerant vr- 
electrons near the Fermi energy which behave as massless Dirac fermion. The valence 
and conduction bands touch conically at two nonequivalent Dirac points, called 
and K- point, which form a time-reversed pair, i.e. opposite chirality. The chirality 
and a Berry phase of vr at the two Dirac points provide an environment for highly 
unconventional and fascinating two-dimensional electronic properties, [2] such as the 
half-integer quantum Hall effect, [3] the absence of backward scattering, [H [5] vr-phase 
shift of the Shubnikov-de Haas oscillations. [6] 

The successive miniaturization of the graphene electronic devices inevitably 
demands the clarification of edge effects on the electronic structures and electronic 
transport properties of nanometer-sized graphene. The presence of edges in graphene 
has strong implications for the low-energy spectrum of the vr-electrons. [3, El E] There 
are two basic shapes of edges, armchair and zigzag which determine the properties of 
graphene ribbons. It was shown that ribbons with zigzag edges (zigzag ribbon) possess 
localized edge states with energies close to the Fermi level. [H [HI [9l [10] These edge states 
correspond to the non-bonding configurations as can be seen by examining the analytic 
solution for semi-infinite graphite with a zigzag edge for which the wave functions of 
the edge states reside on one sublattice only. [7] In contrast, edge states are completely 
absent for ribbons with armchair edges. Recent experiments support the evidence of edge 
localized states. [HI [12] Also, graphene nanoribbons can experimentally be produced by 
using lithography techniques and chemical techniques. [131 [TU [15], [161 [E] 

The electronic transport through graphene nanoribbons shows a number of 
intriguing phenomena such as zero- conductance Fano resonances. |18t IT9] valley 
filtering. [20]. half-metallic conduct ion [2lj, spin Hall effect, [22], and perfectly conducting 
channel. [23] Recent studies also clarify the unconventional transport through graphene 
junctions, quantum point contact and heterojunctions.[l9l [Ml [SS, [261 [23, [SHI [29l [301 
[3ll[32l[33l[3a[35l|33IMll39lll0lliIlll3[ia It is also expected that the edge 
states play an important role for the magnetic properties in nanometer-sized graphite 
systems, because of their relatively large contribution to the density of states at the 
Fermi energy. [Tl[9l[45l[46l[47l[48l[49l[50l[5ll[5a[53[^ Recent studies explore 

the robustness of edge states to the size and geometries [56l [53, |5H], and various edge 
structures and modification. [561 ES] 

Since the graphene nanoribbons can be viewed as a new class of quantum wires, 
one might expect that random impurities inevitably cause Anderson localization, i.e. 
conductance decays exponentially with increasing system length L and eventually 
vanishes in the limit of L oo. However, it was shown that carbon nanotubes 
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with long-ranged impurities possess a perfectly conducting channel. [SI] Recent studies 
show that perfectly conducting channels can be stabilized in two standard universality 
classes. One is the symplectic universality class with an odd number of conducting 
channels. |6H [62l [63] and the other is the unitary universality class with the imbalance 
between the numbers of conducting channels in two propagating directions. [231 EH 165] 
The symplectic class consists of systems having time-reversal symmetry without spin- 
rotation invariance, while the unitary class is characterized by the absence of time- 
reversal symmetry. [66] 

In this paper, we will give a brief overview on the electronic transport properties 
of disordered graphene nanoribbons. In zigzag nanoribbons, the edge states play an 
important role, since they appear as special modes with partially flat bands and lead 
under certain conditions to chiral modes separately in the two valleys. There is one such 
mode of opposite orientation in each of the two valleys of propagating modes, which are 
well separated in /c-space. The key result of this study is that for disorder without 
inter-valley scattering a single perfectly conducting channel emerges introduced by the 
presence of these chiral modes. This effect disappears as soon as inter-valley scattering 
is possible. On the other hand, the low-energy spectrum of graphene nanoribbons with 
armchair edges (armchair nanoribbons) is described as the superposition of two non- 
equivalent Dirac points of graphene. In spite of the lack of well-separated two valley 
structures, the single-channel transport subjected to long-ranged impurities is nearly 
perfectly conducting, where the backward scattering matrix elements in the lowest order 
vanish as a manifestation of internal phase structures of the wavefunction. [63 For multi- 
channel energy regime, however, the conventional exponential decay of the averaged 
conductance occurs. Symmetry considerations lead to the classification of disordered 
zigzag ribbons into the unitary class for long-range impurities, and the orthogonal class 
for short-range impurities. Since the inter- valley scattering is not completely absent, 
armchair nanoribbons can be classified into orthogonal universality class irrespective of 
the range of impurities. 

2. Electronic states of graphene and nanoribbons 
2.1. Tight-binding model and edge states 

There are two typical shapes of a graphene edge, called armchair and zigzag. The two 
edges have 30 degrees difference in their cutting direction. Here we briefly discuss 
the way that the graphene edges drastically change the n electronic structures. [7]. 
Especially, a zigzag edge provides the localized edge state, while an armchair edge does 
not show such localized states. 

A simple and useful model to study the edge and size effect is one of the graphene 
ribbon models as shown in Figs.[l](a) and (b). We define the width of graphene ribbons as 
A^, where N stands for the number of the dimer (two carbon sites) lines for the armchair 
ribbon and by the number of the zigzag lines for the zigzag ribbon, respectively. It is 
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Figure 1. Structure of graphene nanoribbon with (a) armchair edges (armchair 
ribbon) and (b) zigzag edges (zigzag ribbon). The lattice constant is a and N defines 
the ribbon width. The circles with dashed line indicate the missing carbon atoms 
for the edge boundary condition of massless Dirac equation. The disordered region 
with randomly distributed impurities lines in the shaded region and has the length 
L(see text in section 3). Randomly distributed circles schematically represent the 
long-ranged impurities. 

assumed that all dangling bonds at graphene edges are terminated by hydrogen atoms, 
and thus give no contribution to the electronic states near the Fermi level. We employ 
a single-orbital tight binding model for the vr electron network. The Hamiltonian is 
written as, 

H=-tY: c\c, + y: v4c^1 (1) 

{hi) » 

where the operator c\ creates an vr-electron on the site i. denotes the summation 

over the nearest neighbor sites, t the transfer integrals between all the nearest neighbor 
sites are set to be unity for simplicity. This is sufficient to show the intrinsic difference in 
the electronic states originating from the topological nature of each system. The value 
of t is considered to be about 2.75eV in a graphene system. The second term in Eq. ([T]) 
represents the impurity potential, Vi = V{ri) is the impurity potential at a position rj. 
The effect of impurity potential on the electronic transport properties will be discussed 
in the next section. 

Prior to the discussion of the vr-electronic states of graphene nanoribbons, we 
shall briefly review the 7r-band structure of a graphene sheet [68] . To diagonalize the 
Hamiltonian for a graphene sheet, we use a basis of two-component spinor, cj^ =(c^k5 
Cfik)) which is the Fourier transform of (c|g^, cjg^). Let ri, T2, T3 be the displacement 
vectors from a B site to its three nearest-neighbor A sites, defined so that z • t-^ X T2 is 
positive (Fig. 2(a)). z is the normal vector to the graphene sheet. In this representation, 
the Hamiltonian is written as if = J^ic Ck-f^kCk and 

3 

ifk = XI (cos(fc • Ti)a^ + sin(fc • Ti)ay) , (2) 

i=l 

where <t = {a^:, d'y, a^) are the Pauli matrices. Then, the energy eigenvalues are 
iexp{k • Ti)\. Since one carbon site has one vr-electron on average, only 
-E,^-band is completely occupied. 
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Figure 2. (a) Graphene sheet in real space, where the black (white) circles 
mean the A(B)-sublattice site, a is the lattice constant. Here ti — (0,a/\/3), 
T2 = (—a/2, — a/2\/3), and T3 = (a/2, — a/2\/3). (b) First Brillouin zone of graphene. 
K+ = ^(^75), K_ = ^(§,0), r = (0,0) (c) The tt band structure and (d) the 
density of states of graphene sheet. The valence and conduction bands make contact 
at the degeneracy point Kj-. 

In FigsE] (b)-(d), the 1st Brillouin Zone (BZ) of graphene lattice, the energy 
dispersion of vr-bands in the 1st Brillouin Zone (BZ) and the corresponding density 
of states are depicted, respectively. Near the T point, both valence and conduction 
bands have the quadratic form of and ky, i.e. Ey. = ±(3 — 3|fcp/4). At the M points, 
the middle points of sides of the hexagonal BZ, the saddle point of energy dispersion 
appears and the density of states diverges logarithmically. Near the K point of the 
corner of hexagonal 1st BZ, the energy dispersion is linear in the magnitude of the wave 
vector, Ek = ±\/3ta\k\/2, where the density of states linearly depends on energy. Here 
a(= \/3|Ti|(i = 1,2,3)) is the lattice constant. The Fermi energy is located at the K 
points and there is no energy gap at these points, since E^ vanishes at these points by 
the hexagonal symmetry. 

The energy band structures of armchair ribbons are shown in Figsl3] (a)-(c), for 
three different ribbon widths, together with the density of states. The wave number 
k is normalized by the length of the primitive translation vector of each graphene 
nanoribbon, and the energy E is scaled by the transfer integral t. The top of the 
valence band and the bottom of the conduction band are located at /c = 0. It should be 
noted that the ribbon width decides whether the system is metallic or semiconducting 
. As shown in Figj3] (b), the system is metallic when = 3M — 1, where M is 
an integer. For the semiconducting ribbons, the direct gap decreases with increasing 
ribbon width and tends to zero in the limit of very large A^. For narrow non-doped 
metallic armchair nanoribbons, the energy gap can acquire due to Peierls instabilities 
toward low-temperatures. [69] which is consistent with the recent density functional 
theory calculation. [561 EO] 

For zigzag ribbons, however, a remarkable feature arises in the band structure, as 
shown in FigsH] (a)-(c). We see that the highest valence band and lowest conduction 
band are always degenerate at = vr. It is found that the degeneracy of the center bands 
at /c = vr does not originate from the intrinsic band structure of graphene sheet. These 
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Figure 3. Energy band structure E{k) and density of states D(E) of armchair ribbons 
of various widths [(a) = 4, (b) 5 and (c) 30 ]. 



(a) (b) (c) 
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Figure 5. Charge density plot for analytic solution of the edge states in a semi- 
infinite graphite, when (a)A: = 7r, (b)87r/9, (c)77r/9 and (d)27r/3. (e) An analytic form 
of the edge state for a semi-infinite graphite sheet with a zigzag edge, emphasized by 
bold lines. Each carbon site is specified by a location index n on the zigzag chain and 
by a chain order index m from the edge. The magnitude of the charge density at each 
site, such as x, y and z, is obtained analytically (see text). The radius of each circle is 
proportional to the charge density on each site, and the drawing is made for k = 7n/9. 

two special center bands get flatter with increasing ribbon width. A pair of partial flat 
bands appears within the region of 27r/3 < |A;| < tt, where the bands sit in the vicinity 
of the Fermi level. 

The electronic state in the partial fiat bands of the zigzag ribbons can be 
understood as the localized state near the zigzag edge via examining the charge density 
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distribution [71 El El IHl Here we show that the puzzle for the emergence of the edge 
state can be solved by considering a semi-infinite graphite sheet with a zigzag edge. 
First to show the analytic form, we depict the distribution of charge density in the flat 
band states for some wave numbers in Fig|5](a)-(d), where the amplitude is proportional 
to the radius. The wave function has non-bonding character, i. e. finite amplitudes only 
on one of the two sublattices which includes the edge sites. It is completely localized 
at the edge site when = vr, and starts to gradually penetrate into the inner sites as k 
deviates from vr reaching the extend state at A; = 27r/3. 

Considering the translational symmetry, we can start constructing the analytic 
solution for the edge state by letting the Bloch components of the linear combination of 
atomic orbitals (LCAO) wavefunction be e''^("'~^\ e'*^", e'*^*^'^"'"^ on successive edge 
sites, where n denotes a site location on the edge. Then the mathematical condition 
necessary for the wave function to be exact for E = is that the total sum of the 
components of the complex wave function over the nearest-neighbor sites should vanish. 
In Figj5](e), the above condition is e'''("+^) + e'^" + a; = 0, e'''" + e'''^""-^) + y = and 
X + y + z = 0. Therefore, the wave function components x, y and z are found to be 
jj^^ik{n+i/2) ^ ^^gi/c(n-i/2)^ ^'le^''", respectively. Here Dk = -2cos(A;/2). We can thus see 
that the charge density is proportional to Df.^"^~^^ at each non-nodal site of the m-th 
zigzag chain from the edge. Then the convergence condition of \Dk\ < 1 is required, 
for otherwise the wave function would diverge in a semi-infinite graphite sheet. This 
convergence condition defines the region 2tt/3 < \k\ < vr where the flat band appears. 



2.2. massless Dirac equation 



We briefly discuss here the relation between massless Dirac spectrum of graphene and 
low-energy electronic states of nanoribbons. The electronic states near the two non- 
equivalent Dirac points {K±) can be described by 4 x 4 Dirac equation, i.e. 



with 



and 



Hk.pFir) = eF{r) 



Hk-r 



( 

'^{k^ + iky) 



\ 








7(^^ - iky) 










likx ~ iky 








7(4. + iky) 



(3) 



(4) 



Fir) 



\ 



B {r) 

B 



\ 



(5) 



Fi 

F^-{r) J 

Here, kx{ky) is wavevector operator, and can be replaced as — > — iV in the absence 
of magnetic field. 7 is a band parameter which satisfies 7 = ^/d>ta/2. F^^{r) and 



(r) are the envelope functions near K± points for A and B sublattice which slowly 
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vary in the length scale of the lattice constant. We can rewrite the above effective mass 
Hamiltonian by using the Pauli matrices r^^y'^ for valley space {K±) as 

Hk.p = 7 [kicr'^ ® r') + kyia^ ® r')] . (6) 

Here, r° is the 2x2 identity matrix. We can easily obtain the linear energy spectrum 
for graphene as 

e = S'y\k\ with s = ±1, (7) 



and the corresponding wavef unctions with the definition of ^k± 



Fk±Ai Fk±b 



are 



ik-r 



(8) 



Here 

±i0fc _ ± iky 

2.2.1. Zigzag nanoribbons The low-energy electronic states for zigzag nanoribbons also 
can be described starting from the Dirac equation. [TOl [71] Since the outermost sites along 
1^* {N^^) zigzag chain are B(A)-sublattice, an imbalance between two sublattices occurs 
at the zigzag edges leading to the boundary conditions 

(t)K±A{.r^) = 0, (l)K±Bir[N+i]) = 0, (10) 

where r[j] stands for the coordinate at z*^ zigzag chain. The energy eigenvalue and 
wavenumber is given by the following relation, 

e = ±iri - k)e''^ , (11) 



where rj = y/k'^ — e. It can be shown that the valley near k = 37r/2a in Fig. 1(b) 
originates from the l^_|_-point, the other valley at k = — 37r/2a from -K^_-point.[TO| [71] 



2.2.2. Armchair nanoribbons The boundary condition of armchair nanoribbons 
projects and states into F point in the first Brillouin Zone as can be seen in 
Fig.[2](b). Thus, the low-energy states for armchair nanoribbons are the superposition 
of and states. The boundary condition for armchair nanoribbons [71] can be 
written as 

[F+ix,y) + Fxix,y)]\,=o,w = 0, (12) 

[F+ix,y)-Fs{x,y)]U=o,w = 0. (13) 

If the ribbon width W satisfies the condition of ly = {3/2){Nw + l)a with N^j = 
0,1,2,..., the system becomes metallic with the linear spectrum. The corresponding 
energy is given by 

en,k,s = S7\/<T^, (14) 
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where k„ = „..^^"-|. , n = 0, ±1, ±2, . . . and s = ±. The n = mode is the lowest hnear 
subband for metalhc armchair ribbons. The energy gap (A^) to first parabohc subband 
of n = 1 is given as 

A, = 47r7/3(iV^ + l)a, (15) 

which is inversely proportional to ribbon width. It should be noted that small energy gap 
can be acquired due to the Peierls distortion for half-filling at low temperatures [691 170] . 
but such effect is not relevant for single-channel transport in the doped energy regime. 



2.3. Edge boundary condition and intervalley scattering 

Now we discuss the relation between the intervalley scattering and edge boundary 
condition. According to Ref. |5j, the impurity potential can be included in the massless 
Dirac equation by adding the following potential term Uimp described as 



imp 



/ UA{r) 


u'Air)* 







-u'sir) 



u A[r) 


UA{r) 








-u B[r) 


Mij(r) / 



(16) 



with 



'^x{r) =Y,g{r - Rx)ux{Rx) , 

Rx 

« x(r) = J2gir-Rx) e-^^'^-wx [Rx] 
Rx 



[17] 

:i8) 



where ux{Rx) is the local potential due to impurities for X = A oi B. Here g{R) with 
the normalization condition of Z]r5'(-R) = 1 is a real function which has an appreciable 
amplitude in the region where |-R| is smaller than a few times of the lattice constant, 
and decays rapidly with increasing \R\. For convenience we distinguish the impurity 
into two types by the range of the impurity potential: one is long-ranged impurities 
(LRI) if the range of impurity potential is much larger than the lattice constant and the 
other is short-ranged impurity (SRI) if the range of impurity is smaller than the lattice 
constant. 

If only the LRI are present, we can approximate ma(^) = usir) = u{r) and 
^'a(^) = u' sir) = u'{r). In the case of carbon nanotubes and zigzag nanoribbons, 
u' xij"), u'x{r) vanishes after the summation over Rx in Eq. ( |T8i) since the phase factor 
Q-i2K Rx strongly oscillates in the x-direction. This means that the two valleys are 
independent and one can only focus on either or K valley. Thus LRIs do not 
induce the intervalley scattering for zigzag nanoribbons. 

However, this cancellation is not complete in an armchair nanoribbon because the 
averaging over the x-direction is restricted to the finite width of W. This means that we 
cannot neglect the contribution from scatterers particularly in the vicinity of the edges 
to u'x{r). This means that intervalley scattering does not vanish even in the case of 
LRI in the armchair nanoribbons. 
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3. Electronic transport properties 

We numerically discuss the electronic transport properties of the disordered graphene 
nanoribbons. In general, electron scattering in a quantum wire is described by the 
scattering matrix. [66] Through the scattering matrix S, the amplitudes of the scattered 
waves O are related to the incident waves /, 



(19) 



Here, r and r' are reflection matrices, t and t' are transmission matrices, L and R denote 
the left and right lead lines. The Landauer-Biittiker formula [73j relates the scattering 
matrix to the conductance of the sample. The electrical conductance is calculated using 
the Landauer-Biittiker formula, 

G{E) = ^Tritt^) = ^g{E). (20) 
nil nn 

Here the transmission matrix t{E) is calculated by means of the recursive Green 
function method. flEl [71] For simplicity, throughout this paper, we evaluate electronic 
conductance in the unit of quantum conductance (e^/vr^), i.e. dimensionless 
conductance g{E). We would hke to mention that recently the edge disorder effect on 
the electronic transport properties of graphene nanoribbons was studied using similar 
approach. [75l EZj 

3.1. One-way excess channel system 

In this subsection, we consider the conductance of zigzag nanoribbons in the clean limit, 
which is simply given by the number of the conducting channel. As can be seen in Fig. 
EU^a), there is always one excess left-going channel in the right valley {K^) within the 
energy window of \E\ < 1. Analogously, there is one excess right-going channel in the 
left valley {K_) within the same energy window. Although the number of right-going 
and left-going channels are balanced as a whole system, if we focus on one of two valleys, 
there is always one excess channel in one direction, i.e. a chiral mode. 

Now let us consider to inject electrons from left to right-side through the sample. 
When the chemical potential is changed from E = 0, the quantization rule of the 
dimensionless conductance {gK+) in the valley of X+ is given as 

9k+ = n, (21) 

where n = 0, 1, 2, ■ ■ -. The quantization rule in the K_-val\eY is 

9K+=n + 1. (22) 

Thus, conductance quantization of the zigzag nanoribbon in the clean limit near E = 
has the following odd-number quantization, i.e. 

9 = gK++gK_=2n + l. (23) 
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Figure 6. (a) Energy dispersion of zigzag ribbon with N ~ 10. The vaUeys in the 
energy dispersion near k = 27r/3a (fc = — 27r/3a) originate from the Dirac K^{K-)- 
point of graphene. The red-filled (blue-unfilled) circles denote the right (left)-moving 
open channel at the energy i?o(dashed horizontal line). In the left(right) valley, 
the degeneracy between right and left moving channels is missing due to one excess 
right(left)-going mode. The time-reversal symmetry under the intra-valley scattering 
is also broken, (b) Schematic figure of scattering geometry at and K— points in 
zigzag nanoribbons, where a single excess right-going mode exists for K— point. But 
a single excess left-going mode exists for point. Here ric = 0, 1, 2, • • •. 



Since we have an excess mode in each valley, the scattering matrix has some 
peculiar features which can be seen when we explicitly write the valley dependence 
in the scattering matrix. By denoting the contribution of the right valley (-K'+) as -|-, 
and of the left valley {K-) as — , the scattering matrix can be rewritten as 



/ ot \ 



t' 



(It\ 



It 



(24) 



Here we should note that the dimension of each column vector is not identical. Let us 
denote the number of the right-going channel in the valley Kj^ or the left-going channel 
in the valley as ric- For example, Uc = 1 aX E = Eq m Figj^a). Figl6](b) shows the 
schematic figure of scattering geometry for Kj^ and K— points. Thus the dimension 
of the column vectors is given as follows: 



and 



dim{lt) 
dim{I^) 

dim{0~t) 
dim{OJ^) 



ric + 1, 



+ 1, 



dim{lt) 
dim{I^) - 

dim{Ot) 
dim{0^) 



+ 1, 



Tic 



1. 



(25) 



(26) 



Subsequently, the reflection matrices have the following matrix structures. 

Tie ric + l 
^ ^ + 1 / r++ 

ric 



Tie 

r+- 
r 



(27) 
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r'= i 7+ 7- ). (28) 

nc + l\ r'_^ r'__ ) ^ ' 

The reflection matrices become non-square when the intervalley scattering is suppressed, 

i.e. the off-diagonal submatrices (r_| , r |_ and so on) are zero. 

When the electrons are injected from the left lead of the sample and the intervalley 
scattering is suppressed, a system with an excess channel is realised in the K.—- 

valley. Thus, for single valley transport, the r and r'_ are x (n^ + 1) and 

[uc + 1) X matrices, respectively, and t and t'_ are [uc + 1) x (n^ + 1) and 

Uc X ric matrices, respectively. Noting the dimensions of r and r'_ , we find that 

r\ r and r'_ r'\ have a single zero eigenvalue. Combining this property with the 

flux conservation relation [S'^S = SS'^ = 1), we arrive at the conclusion that t t! 

has an eigenvalue equal to unity, which indicates the presence of a perfectly conducting 

channel (PCC) only in the right-moving channels. Note that t' t'l does not have 

such an anomalous eigenvalue. If the set of eigenvalues for t' is expressed as 

{Ti, T2, • ■ ■ , T„^}, that for t tl is expressed as {Ti, T2, ■ ■ ■ , T„^, 1}, i.e. a PCC. Thus, 

the dimensionless conductance g for the right-moving channels is given as 

rtc+l ric 

gK_= J2T^ = l + T.T^, (29) 

i=l i=l 

while that for the left-moving channels is 

9'K_=Y.T^. (30) 

We see that gK_ = g'K_ + 1- Since the overall time reversal symmetry (TRS) of the 
system guarantees the following relation: 

9'k^=9k_, 
9k_ = 9k+, 

the conductance g = gK+ + gK_ (right-moving) and g' = g'j^^ + g'j^_ (left-moving) are 
equivalent. If the probability distribution of {Tj} is obtained as a function L, we can 
describe the statistical properties of g as well as g' . The evolution of the distribution 
function with increasing L is described by the DMPK(Dorokhov-Mello-Pereyra-Kumar) 
equation for transmission eigenvalues. [65] 

In the following, the presence of a perfectly conducting channel in disordered 
graphene nanoribbons will be demonstrated with the help of numerical calculation. 
Recently Hirose et. al. pointed out that the Chalker-Coddington model which possesses 
non-square reflection matrices with unitary symmetry gives rise to a perfectly conducting 
channel. [6l] However, systems with an excess channel in one direction has been believed 
difficult to realize. Therefore disordered graphene zigzag nanoribbons with LRI might 
constitute the first realistic example. It is possible to extend the discussion to generic 
multiple-excess channel model, where the m-PCCs (m = 2,3,---) appear. [65] Such 
systems can be realized by stacking zigzag nanographene ribbons. [78] The electronic 



Electronic transport properties of graphene nanoribbons 



13 



transport due to PCC resembles to the electronic transport due to a chiral mode in 
quantum Hall system. However, it should be noted that the PCC due to edge states in 
zigzag ribbons occurs even without the magnetic field. [791 EQ] 



3.2. Model of impurity potential 

As shown in Fig.l, the impurities are randomly distributed with a density riimp in the 
nanoribbons. In our model we assume that the each impurity potential has a Gaussian 
form of a range d 

V{r,)= "^^Pf-^^"^) (32) 

ro (random) \ ' 

where the strength u is uniformly distributed within the range |m| < um- Here um 
satisfies the normalization condition: 

(jull space) 

Um E exp (-rfd^) /(v/3/2) = no. (33) 

ri 

In this work, we set nimp. = 0.1, uq = 1.0 and d/a = 1.5 for LRI and d/a = 0.1 for SRI. 



3.3. Perfectly conducting channel: absence of Anderson localization 

We focus first on the case of LRI using a potential with d/a = 1.5 which is already 
sufficient to avoid inter-valley scattering. FigI7](a) shows the averaged dimensionless 
conductance as a function of L for different incident energies (Fermi energies), averaging 
over an ensemble of 40000 samples with different impurity configurations for ribbons 
of width = 10. The potential strength and impurity density are chosen to be 
Uq = 1.0 and Uimp. = 0.1, respectively. As a typical localization effect we observe 
that {g) gradually decreases with growing length L (FigJTj). However, {g) converges 
to {g) = 1 for LRIs (FigJTl^a)), indicating the presence of a single perfectly conducting 
channel. It can be seen that (g) (L) has an exponential behavior as 

{g)-l^ exp(-L/0 (34) 

with ^ as the localization length. 

We performed a number of tests to confirm the presence of this perfectly conducting 
channel. First of all, it exists up to L = 3000a for various ribbon widths up to A^ = 40 
for the potential range {d/a = 1.5). Moreover the perfectly conducting channel remains 
for LRI with d/a = 2.0, 4.0, 6.0, 8.0, and Uq = 1.0, Uimp. = 0.1 and A^ = 10. As the effect 
is connected with the subtle feature of an excess mode in the band structure, it is natural 
that the result can only be valid for sufficiently weak potentials. For potential strengths 
comparable to the energy scale of the band structure, e.g. the energy difference between 
the transverse modes, the result should be qualitatively altered. [8lj Deviations from 
the limit (g) —>■ 1 also occur, if the incident energy lies at a value close to the change 
between g = 2n — l and g = 2n + l for the ribbon without disorder. This is for example 
visible in above calculations for E = 0.4 where the limiting value (g) < 1 (FiglTl^a)). 
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Figure 7. L-dependence of the averaged dimensionless conductance, (g) for zigzag 
nanoribbon with N = 10, (a) d/a = 1.5 (no inter-valley scattering), (b) d/a — 0.1 
(inter- valley scattering). Here uq = 1.0, and riimp. =0.1. More than 9000 samples 
with different impurity configuration are included in the ensemble average. 

Turning to the case of SRI the inter-valley scattering becomes sizable enough to 
ensure TRS, such that the perfect transport supported by the effective chiral mode 
in a single valley ceases to exist. In FiglTl^b), the nanoribbon length dependence 
of the averaged conductance for SRIs is shown. Since SRI causes the inter-valley 
scattering for any incident energy, the electrons tend to be localized and the averaged 
conductance decays exponentially, (g) ~ exp(— L/^), without developing a perfect 
conduction channel. 

In this subsection, we have completely neglected the effect of electron-electron 
interaction, which may acquire the energy gap for non-doped zigzag nanoribbon at 
very low-temperatures accompanying with the edge spin polarization. [TJ HSl [21] In such 
situation, small transport gap will appear near E = 0. Since the edge states have less 
Fermi instability for doped regime, the spin polarized states might be less important for 
doped system. 

3.4- Nearly pefect single channel transport in disordered armchair nanoribbons 

Now we turn to the discussion of the electronic transport properties of disordered 
metallic armchair nanoribbons. Figure [H](a) shows the averaged conductance (g) as a 
function of the ribbon length L in the presence of LRI for several different Fermi energies 
E. As we can clearly see, the averaged conductance subjected to LRI in the single- 
channel transport {E = 0.1,0.2 and 0.3) is nearly equal to one even in the long wire 
regime. This result is contrary to our expectation that electrons are scattered even by 
LRI, since wavefunctions at and K- points are mixed in armchair nanoribbons as we 
have already seen in Sec J2.2.2l For multi-channel transport {E > 0.4), the conductance 
shows a conventional decay. The robustness of single- channel transport can be clearly 
viewed from the Fermi energy dependence of conductance for several different ribbon 
lengths L as shown in Fig.[8](b). It should be noted that the energy dependence in the 
vicinity of £^ = is quite different from that in zigzag nanoribbons. The conductance 
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Figure 8. (a) Average conductance (g) as a function of the ribbon lengtli L in 
the presence of long-ranged impurities (LRI) for several different Fermi energies 
E. Conductance is almost unaffected by impurities for single-channel transport 
{E = 0.1,0.2 and 0.3) while it shows a conventional exponential decay for multi- 
channel transport {E > 0.4). Here, N = 14, nimp. = 0.1 and d/a = 1.5. Ensemble 
average is taken over 10''' samples, (b) The Fermi energy dependence of (g) for LRI. 
(c) The same as (a) for short-ranged impurities (SRI). Here, N = 14, ni,„p. = 0.1 and 
d/a = 0.1. 



decays rapidly due to the finite ribbon width effect in zigzag ribbons [23] while the 
conductance around E = remains unity in armchair ribbons (Fig.[8](b)). 

Now let us see the effect of short-ranged impurities (SRI). Figure [Hl^c) shows the 
average conductance (g) as a function of the ribbon length L in the presence of SRI for 
several different Fermi energies E. In this case, the conductance decays exponentially 
even for single-channel transport. This result is similar to that previously obtained in 
zigzag nanoribbons. However, the rate of decay in the low-energy single-channel regime 
{E = 0.1 and 0.2) is slower than that for multi-channel transport regime {E > 0.4) 
in this case. Similar results are obtained in Ref. [75], but in which only short-ranged 
disorder at the edge of ribbon is considered. 



3.5. T -matrix analysis 

The absence of localization in the single- channel region can be understood from the Dirac 
equation including the impurity potential term U\-aip with armchair edge boundary. To 
consider the amplitude of backward scattering, we introduce the T-matrix defined as 

1 

'imp 



T = U- -I- U- 



E-Hn 



imp ~l~ 



(35) 



We can evaluate the matrix elements of f/imp for the eigenstate |n. A;, s) with the 
eigenenergy of Eq. (1141) which can be written as 

/ / s \ \ 



|n, /c, s) 



1 



— S 

-W(n,k) 



(36) 
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with the phase factor 



^-win,k) ^ ^» . (37) 



Here it should be noted that the phase structure in Eq. (!36l) is different between 
and K_ states, and this internal phase structures are critical for the scattering matrix 
elements of armchair nanoribbons as we discuss in the following. Using the above 
expression, we can obtain the scattering matrix element 

(n, k, S|f/i^p|n', k', S') = (SS' + emn,k)-ein\k'))^^ y ^. ^ (gg) 

with 

\ j-W 

V (n, k- n', k') = —— / dx dy e'^^'^"^ )^ 
AW L Jo Jo 

X 'u{r) (e-^(^»-^n'> + c.c.) - (u'{r)e''^''"+''"'^^ + c.c.)] . (39) 

It should be emphasized that Eq. (1381) has the same form as that obtained for carbon 
nanotubes without inter- valley scattering {u'x{r) = 0) [5J. Interestingly, in spite of the 
fact that armchair nanoribbons inevitably suffer from the inter-valley scattering due to 
the armchair edges {u'x{r) 7^ 0), we can express the matrix element for the backward 
scattering as Eq. (l38l) by including u'x{r) into V{n, k; n', k') in Eq. (1391) . This is due to 
the different phase structure between Kj^ and K- in Eq. (l36i) . 

We focus on the single-channel regime where only the lowest subband with n = 
crosses the Fermi level. From Eq. (1381) . the scattering amplitude from the propagating 
state |0,A;, s) to its backward state |0, — A;, s) in the single-channel mode becomes 
identically zero, i.e. 

(0,-A;,s|?7in,p|0,A;,s) = 0. (40) 

Thus, since the lowest backward scattering matrix element of T-matrix vanishes, the 
decay of {g) in the single channel energy regime is extremely slow as a function of the 
ribbon length as we have seen in Fig. [HI However, the back-scattering amplitude in the 
second and much higher order does not vanish. Hence the single-channel conduction is 
not exactly perfect like carbon nanotubes [5], but near/y perfect in armchair nanoribbons. 



4. Universality class 

According to random matrix theory, ordinary disordered quantum wires are classified 
into the standard universality classes, orthogonal, unitary and symplectic. The 
universality classes describe transport properties which are independent of the 
microscopic details of disordered wires. These classes can be specified by time-reversal 
and spin rotation symmetry. The orthogonal class consists of systems having both time- 
reversal and spin-rotation symmetries, while the unitary class is characterized by the 
absence of time-reversal symmetry. The systems having time-reversal symmetry without 
spin-rotation symmetry belong to the symplectic class. These universality classes have 
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been believed to inevitably cause the Anderson localization although typical behaviors 
are different from class to class. 

In the graphene system, the presence or absence of the intervalley scattering affect 
the time reversal symmetry of the system. If the inter- valley scattering is absent, i.e. 
u'x{r) = 0, the Hamiltonian Hq + Uimp becomes invariant under the transformation 
of 5 = — i (cr^ ® T^) C, where C is the complex-conjugate operator. This operation 
corresponds to the special time-reversal operation for pseudospins within each valley, 
and supports that the system has the symplectic symmetry. However, in the presence 
of inter-valley scattering due to SRI, the invariance under S is broken. In this case, the 
time reversal symmetry across two valleys described by the operator T = (cr^ Cg> r^') C 
becomes relevant, which indicates orthogonal universality class. Thus as noted in 
Ref . [72] , graphene with LRI belongs to symplectic symmetry, but that with SRI belongs 
to orthogonal symmetry. 

However, in the zigzag nanoribbons, the boundary conditions which treat the two 
sublattices asymmetrically leading to edge states give rise to a single special mode in each 
valley. Considering now one of the two valleys separately, say the one around k = k^, we 
see that the pseudo TRS is violated in the sense that we find one more left-moving than 
right-moving mode. Thus, as long as disorder promotes only intra-valley scattering, 
the system has no time-reversal symmetry. On the other hand, if disorder yields inter- 
valley scattering, the pseudo TRS disappears but the ordinary TRS is relevant making 
a complete set of pairs of time-reversed modes across the two valleys. Thus we expect 
to see qualitative differences in the properties if the range of the impurity potentials is 
changed. 

The presence of one perfectly conducting channel has been recently found in 
disordered metallic carbon nanotubes with LRI. [61] The PCC in this system originates 
from the skew-symmetry of the reflection matrix, V = — 'r,[6T] which is special to the 
symplectic symmetry with odd number of channels. The electronic transport properties 
such system have been studied on the basis of the random matrix theory. [62l [63] On the 
other hand, zigzag ribbons without inter-valley scattering are not in the symplectic class, 
since they break TRS in a special way. The decisive feature for a perfectly conducting 
channel is the presence of one excess mode in each valley as discussed in the previous 
section. 

In view of this classification we find that the universality class of the disordered 
zigzag nanoribbon with long-ranged impurity potential (no inter-valley scattering) is 
the unitary class (no TRS). On the other hand, for short-range impurity potentials with 
inter- valley scattering the disordered ribbon belongs to the orthogonal class (with overall 
TRS). Consequently we can observe a crossover between two universality classes when 
we change the impurity range continuously. 

However, in the disordered armchair nanoribbons, the special time-reversal 
symmetry within each valley is broken even in the case of LRI. This is because u'x{r) ^ 
as we have seen in Section [231 Thus, irrespective of the range of impurities, the armchair 
nanoribbons are classified into orthogonal universality class. Since the disordered zigzag 
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Figure 9. Summary concerning the universality crossover. With increasing the range 
of the impurity potential, graphene is known to be the orthogonal for SRIs and the 
symplectic class for LRIs. However, zigzag nanoribbons are unitary class for SRIs. 
Armchair ribbons are classified into orthogonal class for all the range of the impurity, 
is the phase coherence length. W is the width of graphene ribbons. 



nanoribbons are classified into unitary class for LRI but orthogonal class for SRI 12^ . 
it should be noted that the universality crossover in nanographene system can occur not 
only due to the range of impurities but also due to the edge boundary conditions. 

The application of a magnetic field enforces the above arguments. In FigJTO| 
the magnetic field dependence of the averaged conductance for disordered zigzag 
nanorribons with (a) LRIs and (b) SRIs is presented. Similarly, the case for armchair 
nanoribbons with (c) LRIs and (d) SRIs. Here we have included the magnetic filed 
perpendicular to the graphite plane which are incorporated via the Peierls phase: 



7i,iexp 



i27r // dl • A , where A is the vector potential. Since the time- 
reversal symmetry within the valleys for zigzag ribbons is already broken, the averaged 
conductance {g) for LRIs (absence of intervalley scattering) is quite insensitive to the 
application of magnetic field as can be seen in Fig JTOT a). This is consistent with the 
behavior in the unitary class. For higher energies, weak magnetic field dependence 
appears due to the intervalley scattering. For all the other cases weak magnetic filed 
improve the conductance, i.e. weak localization behavior which is typical for the 
orthogonal class. 



5. Summary 

In this paper, we have presented a brief overview on the electronic and transport 
properties of graphene nanoribbons focusing on the effect of edge shapes and impurity 
scattering. Concerning transport properties for disordered systems the most important 
consequence is the presence of a perfectly conducting channel in zigzag nanoribbons, i. e. 
the absence of Anderson localization which is believed to inevitably occur in the one- 
dimensional electron system. The origin of this effect lies in the single-valley transport 
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Figure 10. The magnetic field dependence of the averaged conductance for disordered 
zigzag nanoribbons with (a) LRIs and (b) SRIs. Similarly, the case for armchair 
nanoribbons with (c) LRIs and (d) SRIs. cj) is the magnetic flux through a hexagon 
ring which is measured in units of ch/e. 



which is dominated by a chiral mode. On the other hand, large momentum transfer 
through impurities with short-range potentials involves both valleys, destroying this 
effect and leading to usual Anderson localization. The obvious relation of the chiral mode 
with time reversal symmetry leads to the classification into the unitary and orthogonal 
class depending on the range of impurity potential. On the other hand, in spite of the 
lack of well-separated two valley structures, the single-channel transport subjected to 
long-ranged impurities shows nearly perfect transmission, where the backward scattering 
matrix elements in the lowest order vanish as a manifestation of internal phase structures 
of the wavefunction. These results are in contrast with the mechanism of perfectly 
conducting channel in disordered zigzag nanoribbons and metallic nanotubes where the 
well separation between two non-equivalent Dirac points is essential to suppress the 
inter-valley scattering. 
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